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A One-step Approach to Computing a Polytopic Robust Positively Invariant 

Set 


Paul Trodden, Member, IEEE 


Abstract 

A procedure and theoretical results are presented for the problem of determining a minimal robust 
positively invariant (RPI) set for a linear discrete-time system subject to unknown, bounded disturbances. 
The procedure computes, via the solving of a single LP, a polytopic RPI set that is minimal with respect 
to the family of RPI sets generated from a hnite number of inequalities with pre-dehned normal vectors. 


Index Terms 

Linear systems; Uncertain systems; Computational methods; Optimization; Invariant sets 


I. Introduction 

We consider the problem of finding, for the discrete-time, linear time-invariant system, 

x'^ = Ax -f w, (1) 

a robust positively invariant (RPI) set. That is, a set 7^ C with the property 

Ax + w E'lZ^'ix eW . (2) 

In this problem, x G is the current state and x+ its successor. The disturbance w E MA 
is unknown but lies in a polytopic (compact and convex) set W that contains the origin in its 
interior. 

Robust or disturbance invariant sets are important in control, and their theory and computation 
have attracted significant attention; see, for example, |[I1-|111 and references therein. One set that 
is of particular interest is the minimal RPI (mRPI) set—that is, the RPI set that is smallest in 
volume among all the RPI sets for a system—which is also the set of states reachable from the 
origin in the presence of a bounded disturbance. This set is an essential ingredient in many robust 
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control algorithms. For example, in tube-based robust model predictive control (MFC) [j5||, an RPI 
set is used to guarantee robust stability and feasibility in the presence of bounded uncertainty; 
moreover, since the constraints in the MFC optimization problem are tightened according to the 
size of the RPI set, then the smallest RPI set {i.e., the mRPI set) is desirable. However, computing 
an exact representation of the mRPI is generally impossible (except for special instances of A, as 
identified later), and instead an approximation is usually sought. A seminal contribution in this 
regard is [|3l, which proposes a method for computing an abitrarily close outer-approximation 
to the mRPI set, which is itself RPI. 

The essence of the problem of computing exactly the mRPI set is that this set is, in general, 
not finitely determined. Methods for computing approximations to the mRPI set, including dS], 
rely on finding finite representations of the set. Recently, in the context of tube-based MFC, (H 
introduced and studied the notion of a polytopic RPI set defined by a finite number, r, of a- 
priori selected linear inequalities. For a non-autonomous system = Ax + Bu + w controlled 
by a continuous positively homogeneous control law, u = k{x), the authors showed that the 
RPI set dynamic condition Q has an equivalent representation as r functional inequalities. It 
was established that a fixed-point solution to the functional equation corresponds to an RPI set 
that is minimal, in volume, with respect to the entire family of RPI sets defined by the pre¬ 
selected inequalities, and is an invariant outer-approximation to the mRPI set. To compute this 
set, the authors of [jH give an iterative procedure, based on solving a sequence of LPs, for which 
convergence is guaranteed. 

In this note, we adopt the notions of [[61 and specialize their results to the case of the 
linear autonomous system ([U) (alternatively, the linear non-autonomous system with linear state 
feedback control law) in order to develop a one-step approach, based on solving a single LP, to 
the computation of the smallest RPI set defined by a pre-selected system of inequalities. Though 
simple, to the author’s knowledge this has not appeared in the literature, although there are 
related results; for example, it is known that checking the invariance of an existing poly tope 
is an LP [|2l. On the other hand, the ability to synthesize a near-minimal RPI set by solving a 
single LP potentially paves the way for robust control methods that re-compute the disturbance 
invariant sets on-line, as done in, for example, the recently developed “plug-and-play” approach 
to distributed MPC Q. 

The proposed approach differs to the one of O in one important assumption: the number 
and normal vectors of the inequalities that represent the RPI set are, as in [H, defined a priori. 
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while in Q both are unknown until termination of the algorithm. This a-priori definition, first 
proposed and studied by [|6]|, has two eonsequenees: firstly, the RPI set obtained is not neeessarily 
the mRPI set, or even an abitrarily elose outer-approximation (as it is in Q); however, it is the 
smallest RPI set that ean be represented by the finite number, r, of ehosen inequalities with 
normal veetors {P^ = To make a elear distinetion, in this note we term this the 

(F,r) -mRPI set when the number of ehosen inequalities is r and the matrix of normal veetors 
(the left-hand side of the defining system of inequalities) is P. Seeondly, the method of jSl 
involves solving a sequenee of LPs and then eomputing a Minkowski summation, but here only 
the solving of a single LP is required. The development of the procedure here comprises two 
steps, the enumeration of which also serves to clarify the contribution of this note with respect 
to 0: first, we show that, for the studied linear autonomous system dU), the fixed-point solution 
to the functional equation, which [Q showed is guaranteed to exist, is in fact unique. Secondly, 
we show that the corresponding RPI set—which [O proved to be minimal with respect to the 
family of RPI sets represented by (P, r) —can be computed via a single linear program (LP), as 
an alternative the iterative sequence of LPs proposed by dH. 

Another method that uses a single LP to compute a disturbance invariant set is the optimized 
robust control invariance approach of [j4|, applicable to the linear non-autonomous system = 
Ax-\-Bu + w. Because a robust control invariant (RCI) set—and the associated control policy—is 
obtained, then this subsumes the robust positive invariance (where di fixed control law is assumed) 
considered here. However, that approach optimizes over only those control policies that guarantee 
a finitely determined set, achieved by employing a relaxed variation of the assumption, for ([T]), 
that C aW for some a G [0,1) and finite integer k. In this note, the assumption that 

A has eigenvalues inside the unit circle is required, which is different to the assumption used 
for finite determination of RCI sets in dH, but weaker than the assumption required for finite 
determination of the mRPI set for ([1]). 

The organization of this note is as follows. First, in Section [III it is shown that for the 
system ([T]), the fixed-point solution is, under suitable assumptions, unique. Subsequently, in 
Section Uni it is shown that the (P, r)-mRPI set for dU) may be computed via a single LP. 
Finally, examples are given in Section |IV] to illustrate the practicality of the proposed approach, 
before conclusions are made in Section 0 

Notation: The sets of non-negative and positive reals are, respectively, Mo+ and M+. For 
a, 6 G M"', a < 6 applies element by element. A matrix M is non-negative, denoted M > 0, if 
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Mij > 0 for all i and j. XX is the scaling of a set A” by A G M, defined as {Ax : x G X}. 
AX denotes the image of a set A C MA under the linear map A : M” i—)■ and is given by 

{Ax : X G A}. The support function of a set A is /i,(A, v) = sup{x^x : x G A}. A polyhedron is 
the convex intersection of a finite number of halfspaces, and a poly tope is a closed and bounded 
(hence compact) polyhedron. 

II. Existence and uniqueness of a (P, r)-MRPI set 
For the system ([U), we consider the case of a poly topic disturbance set 

W = {w : Fw < g}, (3) 

where F G RPX", g G and make the following two standing assumptions. 

Assumption 1: The set W contains the origin in its interior. 

Assumption 2: The eigenvalues of A are strictly within the unit circle. 

The former assumption requires that g G M+. The latter assumption implies, as shown in [[T]|, 
that for a given compact disturbance set W there exists a compact RPI set, 7Z, for the system ([U), 
satisfying ([21). 

Assumption 3: The RPI set P is a polytope that contains the origin in its interior. 

Note that Assumptions [U and [3] imply that the support functions to W and IZ, respectively, 
are positive—a key technical property that will be used in this note to establish the existence 
and uniqueness of the RPI set that we aim to compute. 

In this note, following [jH, we consider the RPI set constructed from a finite number, r, of 
inequalities with pre-defined normal vectors. That is, TZ = TZ{q), defined as 

n{q)^ [zEMA ■. Pz<q], (4) 

where P E |Pj^ : i E (1,..., r}} spans R", P, is the Ah row of matrix P, and q E Rg.,.. 

The left-hand side of the inequalities—the matrix P—is to be chosen a priori by the designer. 
The following result, which is an application of Farkas’ Femma, establishes basic conditions on 
the matrices A, P and P for the existence of an RPI set for the system ([T]) given the disturbance 
poly tope dU). 
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Theorem 1 (Adapted from Hennet and Castellan Suppose Assumptions hold. Then 
the set Tl(q) with some q = q h robust positive invariant for the system ([T]) if and only if there 
exist non-negative matriees H G and M G sueh that 


HP = PA 

(5a) 

MF = P 

(5b) 

Hq + Mg < q 

(5c) 


We will assume that P is ehosen so that an RPI set exists: 

Assumption 4: For the ehosen P, and the system (A, W), there exists a g G sueh that (|2]) 
holds for all x G Tl{q). 

Remark 1: While Assumption |4] may appear strong, it is needed to narrow the elass of matriees 
that we consider to those that admit an RPI set. However, the procedure presented in the next 
section includes a easy certification of existence of an RPI set for a chosen P: if an RPI set exists, 
the (P, r)-mRPI set is returned. If no RPI set exists, the optimization problem is unbounded. 

The authors of [[6|| show—in the more general setting of a linear non-autonomous system 
controlled by a positively homogeneous state-feedback control law—that RPI condition dB) is 
equivalent to the functional inequality 

c{q) + d<b{q), (6) 

where, for i = bi{q) = h{TZ{q), P(^), Ci(q) = h(A'R(q), P^), di = h(W,P^). That is, 

the set inclusion requirement is replaced by support function inequalities, which is a standard 
technique dUl. Note that b(q) may be different to g; for example, in the case of redundant inequal¬ 
ities defining P(g). The topological properties of these functions described in the following two 
lemmas are essential to establishing existence and uniqueness of the fixed-point solution to ©• 

Lemma 1 (Adapted from Proposition 1 of ^): Suppose that Assumptions [I]-[3] hold. Then the 
functions b: Mg_,_ i-G- Mq_,_, c: Mq.,. eg- Mg_,_ are continuous and monotonically non-decreasing; that 
is, b(ai) < b{a 2 ) for oi < 02 . Also, d G M^. 

Lemma 2 (Positive homegeneity ofb, c): Suppose Assumptions [2] and [3] hold. Then the functions 
b{-) and c(-) are positively homogeneous; that is b(Xa) = Xb{a) for A > 0, with a similar 
expression for c(-). 

Proof: Consider bfXa) = h{Ll{Xa),Pi) for some a G Mq.,., A > 0 and i G 
By definition of TZf), Tl{Xa) = XTl(a). Thus, h{7l{Xa),Pi) = h(X7l{a),Pi) = Xh {71(a), Pf), 
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for A > 0, where the latter equality follows directly from the definition of the support func¬ 
tion im. Hence, bi{Xa) = Xbi{a), therefore b{Xa) = Xb{a). Positive homogeneity of c(-) may be 
established by the same arguments. ■ 

The next result, which concerns the existence of a fixed-point solution to ®, was established 
by [HI in the setting of a linear non-autonomous system controlled by positively homogeneous 
state-feedback control law, and hence immediately applies to the more specialized case considered 
in this note. 

Theorem 2 (Theorem 1 of ^): Suppose Assumptions [I]-[3] hold. Let Q = {g G Mq_,_ : 0 < 
q < q}- Then, (i) for all q G Q, c{q) + d ^ Q and (ii) there exists at least one q* & Q satisfying 
c{q*) + d = b(q*) = q* if and only if Assumption 01 holds. 

Remark 2: The necessity and sufficiency of Assumption 0] follows by definition. In particular, 
if Assumption 0] does not hold, then there does not exist an RPI set for the system (A, W) with 
the chosen P. 

Remark 3: Note that, in view of the assumptions on g and the properties of 6(-), c(-), and d, 
a fixed-point solution q* must be strictly positive. 

With respect to computing a fixed-point solution, the sequence generated by the iterative 
procedure = c{q^^) -f d, with = 0, converges to the fixed-point solution q* with the 
smallest 1-norm value, ||g*||i [HI Theorem 2]. As the following result states, the corresponding 
set 7l{q*) is RPI, and, in fact, is the minimal (smallest volume) RPI set over the family of RPI 
sets defined by the r inequalities with left-hand side P. 

Lemma 3 (Corollary 1 of Tl(q*) = where 

S = {Tl(q) : q G H} , and "H = {g G Mo+ • c(q') + d< b{q)] 

For convenience, we define this set Tl(q*) as the (P, r)-mRPI set. 

Definition 1 ({P,r)-mRPI set): The (P, r)-mRPI set for system ([B is TZ{q*) where q* = 
b(q*) = c(q*) + d. 

In this note, we propose an alernative to the iterative procedure of [HI- To this end, the next 
result shows that the fixed-point solution to ® is, in fact, unique. This result is then exploited 
in Section [Till wherein the problem of finding the fixed-point solution is cast as an LP. 

Theorem 3 (Uniqueness of fixed-point solution): Suppose Assumptions [I]-0| hold. Then there 
exists a unique q* G satisfying c{q*) + d = b{q*) = q*. 
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Proof: Existence is established by Theorem [2l so it remains to show that q* is unique. Let 
l{q) = c{q)+d — h{q) and f{q) = h{q)—q. Finding the fixed-point solution c{q*) + d = b{q*) = q* 
is equivalent to finding q* such that l{q*) = f{q*) = 0. Suppose there exist q^ G and G 

such that Z(g^) = /(g^) = 0, Z(g^) = /(g^) = 0, and g^ 7 ^ g^, i.e., q^ — q^ f 0. There are two 
possibilities: 

(i) qf>q] for at least one i G {1, ..., r}, with g| < q] otherwise; 

(ii) g^ < g^, with qf < q} for at least one i G {1,.. ., r}. 


Consider case (i) Let 


a = mm 
2=1...r 


qt 


qi 


^ ^ ^ > 0 


ql 


qp 


Strict positivity follows from the discussion in Remark [3l Since qf > q} for at least one i, 
then a < 1. Let s = aq^ < q^. It follows, from positive homogeneity of b{-) and the fact that 
6 (g^) — g^ = 0, that /(s) = b{s) — s = b{aq^) — aq^ = a{b{q^) — g^)) = 0. Similarly, 

l{s) = c(s) + d — b{s) 

= c{aq^) + d — b{aq^) 

= ac{q^) + d — ab{q^) 

= a(c(q^) -b(q^)) +d 
> 0 

where the second line follows from the positive homegeneity of c(-) and b(-), and the strict 
inquality with zero follows from c(g^) — b(q^) = —d, a < 1 and d > 0. Now, by definition of a, 
and since a < 1, then s < g^ with Sp = q^. For the same p, we have /p(g^) = bp{q^) — g^ = 0, 
fp{s) = bp{s) — Sp = 0, and, since s < g\ then bp{s) < bp{q^). In fact, bp{s) = bp{q^), as we 
have already shown that bp{s) = Sp = q^. We also have /p(g^) = Cp{q^) + dp — bp{q^) and lp{s) = 
Cp{s) + dp — bp{s). Because bp{q^) = bp{s) and Cp{s) < Cp(g^), it follows that lp{s) < lp{q^). But 
then 0 = lp{q^) > lp{s) > 0 , and we have a contradiction: therefore, we conclude that case [(I) 


cannot hold, and either case |(ii)| holds or g^ = gk Now consider case |(ii)[ and its equivalent 
statement: g^?^ > qf for at least one i G with g] > g| otherwise. Following the same 

set of arguments, starting with the opposite definitions of a = minimi^^{qf /q}} and s = aq^, 
we find that that case |(ii)| cannot hold either. Therefore, q^ = q^ = q*, and the solution is unique. 
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III. Computing the (P, r)-MRPI set via a single LP 

The problem of computing the (P, r)-mRPI set is that of finding the q that satisfies the 
functional inequality (RPI condition) (6) while attaining the smallest value of Uglli. The results 
in the previous section show that this q in fact satisfies (6) with equality; it is the fixed-point 
solution q*. Therefore, the problem of finding q* is 

q* = argmining'll^ : c{q) + d< b{q)} (7) 

Q 

This is not tractable, as, by the definitions of b{-) and c(-), the constraints are maximization 
problems involving the optimization variable: 

max {PiAx : Px < q} + max {PiW : Fw < g} 

< max {PiX : Px < q} 

for i = 1... r. However, by noting that the fixed-point solution is unique, we may replace the 
problem of (|7]) with the maximization problem 

q* = argmax{||g||^ : c{q) + d = b{q)} 
q 

This problem then easily converts to a linear program, as shown by the following. Introduce 
auxiliary variables G M” and ce® G MA for each RPI inequality i G Then, noting 

that q = b{q) = c{q) + d at the desired fixed-point solution, eliminate q and b{q), leading to the 
problem 

r 

F: q* = c* + d*, where (c*, d*) = arg max Cj -f di (8) 

subject to, for alH G {1,..., r}. 


Ci < PiAC, 

(9a) 

PC <c + d, 

(9b) 

di < Piu\ 

(9c) 

Fu'' < g. 

(9d) 


In this problem, maximizing each c* subject to constraints (l9al) and (l9bl) represents finding the 
vector of support functions to ATZ. Constraint (l9bl) represents Px < b{q), with the condition 
c{q) + d = b{q) enforced. Constraints (l9^ and (l9dl) represent finding d, the vector of support 
functions to W. 
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Remark 4: Note that, by definition, di = h(W, is constant and does not depend on 
q. Therefore, d could be computed prior to solving P, by solving a sequence of LPs, before 
entering the optimization as a parameter. However, our aim is to formulate a single LP (a one- 
step procedure) that computes, simultaneously, d, c and hence q. 

Note that each di and Ui is bounded, via (|9^ and (l9dl) and the assumptions on W. Further note 
that this problem always has a feasible solution, since one can choose, for example, q = d* = 0 
and = cu* = 0. The question, then, is whether an optimal solution exists, or the problem is 
unbounded. To this end, we require the following result, which specializes Theorem [T] to the 
fixed-point solution. 

Proposition 1: Suppose Assumptions [BBl hold. A vector q* satisfies the fixed-point relation 
c{q*) + d = b{q*) = q* if and only if there exist non-negative matrices H G and M e R^^^ 
such that 


PA 

(lOa) 

P 

(10b) 

* 

q 

(10c) 


Proof: Consider the ith element of each of c(g*), d and b{q*), defined by the (primal) LPs 


'*) = maxjPjAx : Px < q*} 

(11a) 

di = max\^PiW : Fw < 

(11b) 

'*) = max|PjX : Px < q*} 

(11c) 


If Assumptions [I]-l4] hold, then by the previous results there exists a q* satisfying the fixed-point 
equation. Moreover, each of the terms in (fTTI) is well defined, which is the case if and only if 
each LP is feasible and attains an finite optimum. Therefore, by weak duality, the dual of each 
LP 


cfq*)-. mm{hj q* : hj P = PiA, hi > 0}, 
di'. g : mjF = Pi,mi > 

bfq*): mm{yjq* : yjP = Pi,yi > O}, 
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is feasible. Examining these dual problems, dual feasible solutions exist if and only if there exist 
non-negative hi ^ W, rrii G W, yi G R'” sueh that 

hJP = P,A, 

mjF = Pi, 

yJP = Pr- 

Applying strong duality, which holds in view of the previous arguments, to each of the three 
LPs 


Ciiq*) = hjq*, 

di = mjg, 

Hq*) = yjq*- 


Collecting all rows i = 1.. .r, 

c{q*) = Hq* 
d = Mg, 
b{q*) = Yq*, 

where HP = PA, MF = P, YP = P. Therefore, it follows that if the fixed-point equation 

c{q*) + d = b{q*) = q* 

is satisfied, then so are the conditions (fT^ : conversely, if (fT^ are satisfied, then so is the 
fixed-point equation. ■ 

Then the main result of this section follows. 

Theorem 4: Suppose Assumptions [I]-[3] hold. If P satisfies Assumption |4l then problem 
P admits an optimal solution corresponding to the fixed-point solution q*. Otherwise, P is 
unbounded above. 

Proof: We use duality to prove the theorem. Our goal is to prove that the optimal solution 
to P satisfies the conditions (fTOl) . for some non-negative H and M, if and only if Assumption |4] 
holds, and that P is otherwise unbounded. Since the primal LP problem P is known to be feasible, 
it suffices to show that the dual problem is feasible—and the solution is as claimed—if and only 
if P satisfies Assumption HI on the other hand, if the dual is infeasible, then by weak duality 
the primal problem P is unbounded. 
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The dual problem is 

r 

D: min 

subjeet to, for alH G {1,..., r}, 

r 

k=l 

k=l 

F^/i* - A^P^Xi = 0, 

Xi, z/i > 0 


( 12 ) 


(13a) 

(13b) 

(13c) 

(13d) 

(13e) 

(13f) 


where e M, /r* G r'* G M, ? 7 * G are the dual variables associated with constraints (l9al)- 
(l9dl) respectively. 

We first suppose the dual problem D is feasible. From (I13al) and (I13bl) . Xi = Ui = l + X]I=i 
for alH = 1,..., r, where /rf is the i\h element of jjA' G M’’. From this and (I13cl) . (ll3dL it follows 
that 


PiA = 


P^ = 


(f 


i\T 


1 + ELi/^f 
ivV 


P, 


F, 


1 +ELi/^r 

where the division is permitted since — 0- Collecting all rows i = 1... r of F, it 

follows that di feasible solution to D satisfies (llOal) and (II Obi) with Hij = (1 + Y7k=i bi) P 0^ 

j = 1.. .r, Mij = T]j/{1 + Yl\=ilA) ^ 0’ i = 1 • ■ - P; therefore, H and M are non-negative 


matrices. 

Now we study the optimal solution to D. Since the primal problem P is known to be feasible, 
and we assumed D to be feasible, then by strong duality (which holds regardless of the feasibility 
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of ©) the optimal solutions to P and D are attained and equal in objeetive value. So, applying 
eomplementary slaekness to (l9al) and (|9^ . 


r 


5^A'(c'-PAr) = o 


i=l 


r 




where * denotes a variable in the optimal solution. Sinee eaeh term in these sums is non-positive, 



for i = 1,..., r. Morever, because (by (I13al) and (I13bl) l A* > 0 and u* > 0, then, at the optimum. 



Hence, 



(14) 


Now consider the inequality (l9dl) . Suppose Fcu^* < g for some i G {1,..., r} (i.e., Fju'^* < gj for 
all j = 1.. .p). Complementary slackness implies that = 0 which in turn implies (from (I13dl) . 
assuming that Pj is not trivially all zeros) that u* = 0; but u* > 1 by (I13bl) . which is a 
contradiction. Hence, there must exist a subset K, C {1,... ,p} of active constraints for which 
FkU^* = gk for fc G /C. But for any j ^ /C, = 0. 

Similarly, consider the inequality (l9bl) . By complementary slackness, if P^'"* < c* + d* then 
jjp* = 0. By (I13cl) . this implies that Pj\* = 0. There are two cases to consider: (i) if any 
elements of Pi A are non-zero then A* = 0; (ii) if Pi A = 0 then A* > 0 is permitted. We leave 
case (ii) for now and consider (i) first. A* = 0 contradicts (I13al) . which requires A* > 1. Hence, 
there must exist a subset J C of active constraints for which = c* + d* for 
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j G J. But for any k ^ J, /i” = 0. As a consequence of the preceding arguments, (fT4l) may 
be re-written as 


c* -f d* = 


1 + ZlLi A 





k&K 


VTF,0J^ 


= Hi[c* -f d*) + MiQ 

where Hi is the Ah row of H and Mj is the Ah row of M. The second line follows because 

Hij = 0 for j ^ JT” and Mik = 0 for /c ^ /C, while = c* + d* for j E J and Ffccu** = Qk 

for k E )C. 

Now case (ii). If PiA = 0 then c* = 0. Moreover, A* > 1 is permitted, so the same 
contradiction is not constructed. Then, however, either < d*, hence fA* = 0, or P,^** = d*, 

with lA* > 0, for j G JT” C {1,..., r}, and /i” = 0 for all k ^ J. Either way. 


c* d* — Hii^c* + d*) + MiQ 


as before. 

Finally, collecting all rows i = 1. .. r, 

H{c* + d*) + Mg = c* + d* 

which is the third condition in (fTOl) . This establishes that the solution to P, if it is attainable, 
satisfies the conditions (fT^ for it to be the fixed-point solution. It is attainable if and only if 
the dual problem D is feasible. Therefore, it remains to show that the D is feasible if and only 
if Assumption in holds. 

First, necessity of Assumption HI Suppose Assumption |4] is not satisfied, but the dual D 
is feasible. By definition, if Assumption H] is not satisfied then for the chosen P and system 
(A,W) there does not exist a q satisfying the functional inequality ®. Therefore, there exists 
no q* satisfying the functional equation and, by Proposition [H the conditions (fTOl) . However, the 
attainable optimal solution to P and D satisfies (ITOl) with non-negative H and M, as has been 
shown. Therefore, we have a contradiction, and conclude the optimal solution is attainable, and 
D is feasible, only if Assumption |4] holds. 

Second, sufficiency of Assumption HI Writing the primal constraints dH) in the form Ax < b, 
where x is the vector of primal decision variables, it follows that the dual constraints (fTSl) may 
be written in the form A^y = c, y > 0, where y is the vector of dual variables and c is the 
coefficients vector in the vectorized form, c^x, of the objective function ([8]). By Farkas’ Femma, 
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a feasible solution to A^y = c, y > 0 exists if and only if Ax > 0 c^x > 0. Hence, 

we aim to show that, if Assumption |4] holds, then for all x satisfying Ax > 0 we also have 
c^x > 0. The system Ax > 0 may be written in terms of the primal variables as 

Ci > PiAC 
PC>c+d 
di > Piu’' 

fJ > 0 

for i = 1. .. r. If Assumption m holds, then HiP = PiA and MiF = Pi for some non-negative 
Hi and Mi. Substituting into the system Ax > 0, 

c, > H,PC 
PC >c+d 
di > MiFoj'^ 

fJ > 0 , 

from which it follows that dj > 0 and q > Hi{c + d), hence c > He. But we also have that, 
if Assumption (jH) holds, then there exists some q G for which 0 < Hq < Hq + Mg < q. 
Applying recursively, 0 < H'^q < Hq < q, H^ > 0 because H > 0, and therefore lim„^oo > 
0, if the limit exists. In fact, because HP = PA, the nullspace of P is A-invariant and P has 
rank n, then the eigenvalues of H are are subset of the eigenvalues of A; hence, lim„_^oo H^ = 0 
because p{A) < 1. Then c > He > limn_ 5 .oo H"-e = 0. Consequently, c^x = Y7i=i + di > 0. 
Therefore, D is feasible if Assumption |4] holds. ■ 


IV. Examples 


We consider the non-autonomous system 


X = 


1 1 


0.5 


X + 


0 1 


1 


u + w, 


(15) 


with te G W = {w G : ||te||oo < 0.1}. This is converted to the linear autonomous system (O 
by use of a state feedback control law u = Kx. 
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A. Computation of {P,r)-mRPI from selected inequalities 

First, we use the feedbaek matrix K = [—0.4345, —1.0285], eorresponding to the infinite- 
horizon LQR solution with cost matrices Q = I and R = 1. Note that in this example the mRPI 
set is not finitely determined, and therefore an approximation is required. 

Figure Ufa)] shows the (P, r)-mRPI sets generated from r = 6, 20 and 48 inequalities, wherein 
the ith row of P is designed as 


Pi 


sin 




COS 




(16) 


i.e., so that Px < 1 is the r-sided regular polygon. Also shown is the outer approximation to 
the mRPI, which is itself RPI, computed using the algorithm of [|3|| and a tolerance e = 10“"^. 
This set, termed the e-mRPI set, is defined by 48 non-redundant inequalities. 

Figure [ JEb)] shows a similar comparison using K = [—0.0796, —0.4068], obtained as the LQR 
solution with Q = I and R = 100. Now the e-mRPI (e = 10“"^) comprises 172 non-redundant 
inequalities, while the (P, r)-mRPI sets computed using the proposed method are shown for 
r = 20, 60 and 172, again using (fT^ for P. 

Table H compares the computation times and number of operations for computing the (P, r)- 
mRPI with those for obtaining the e-outer approximation using the algorithm of |l3l . For the latter, 
the Multi-Parametric Toolbox v3.0 [fTOll was used for set operations, with CPLEX 12.6 as the 
LP solver for support function calculations. For the (P, r)-mRPI set computations {i.e., solving 
the LP), CPLEX 12.6 was used as the LP solver. The platform was a 64-bit Intel Core 17-2600 
at 3.40 GHz with 8 GB RAM. Times are reported as the mean elapsed time over 100 runs. 

Comparison was also made with the iterative procedure of [0 for computing the (P, r)-mRPL 
The iterative procedure is 

gfc+i = c(gfc) -f d with go = 0 


for which g^ — )■ g* as fc — )■ oo. This was implemented in Matlab using the MPT v3.0 ifTOl for 
support function calculations (with CPLEX 12.6 as underlying EP solver). The function c(-) was 
evaluated element by element at each iteration; that is, as r separate support function calculations. 
Eor the simplest case considered of K = [—0.4345, —1.0285] and r = 6 (the first row of TableU), 
the number of iterations to convergence (of Ig^+i — g^l to within a chosen tolerance of 10“®) was 
34, which included the solving of 238 EPs and took a mean total time of 1.7 seconds. At the 
other end of the scale, for the most difficult problem considered (K = [—0.0796, —0.4068] and 
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(a) K = [-0.4345,-1.0285] 



Xl 


(b) K = [-0.0796, -0.4068] 

Fig. 1. Comparison of (P, r)-mRPI and e-mRPI sets for the system l ll5t with different feedback matrices K. 


r = 172), the iterative procedure required 70 iterations, the solving of over 12000 LPs, and took, 
on average, 90 seconds. While these times can, of course, be shortened by using optimized code, 
the intention here is merely to report the times obtained using standard computational tools. 

B. Re-computing the {P,r)-mRPI set given P 

An interesting use of the method is when an RPI set for the system is available, but is desired 
to be re-computed or modified; for example, if the disturbance set changes. Potential applications 
of this include “plug-and-play” tube-based approaches to distributed MPC, wherein a dynamic 
subsystems’ disturbance set evolves over time as other subsystems are added to and removed 
from the system of coupled subsystems [|71; in such situations, one needs a new RPI set that 
takes into account the latest disturbance set. One could re-compute from scratch a new RPI set. 
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but it may be advantageous, in the interests of computation time, to modify an existing RPI set 
instead. In the context of the approach proposed here, the P matrix of the known RPI set may 
be used as a basis for computing the new RPI set. 

For the system (fT5l) with K = [—0.4345, —1.0285] and W = {tc G : |tc|oo < 0.1}, the P 
matrix is obtained as that of the e-mRPI set. For e = 10“^, this comprises 48 inequalities. Now 
suppose the disturbance set enlarges to 


W 


w G 


-0.3 


0.1 

1 


<w < 



-0.4 


0.2 

I 


Figure [2] shows the (P, r)-mRPI and e-mRPI sets based on the new disturbance set, using 
for the former the P matrix from the old e-mRPI set. The (P, r)-mRPI set, computed in 0.03 s 
using the proposed method, is visually indistinguishable from the new e-mRPI set. 


V. Conclusions 

A procedure for computing a polytopic robust positively invariant set for a linear uncertain 
system has been presented. The method, which requires the solution of a single LP, obtains the 
an RPI set that is the smallest among those represented by a finite number inequalities with 
pre-defined normal vectors, and offers an alternative method of computation to the iterative 
procedure of [0. Existence and uniqueness of a solution has been established. The practicality 
of the approach has been demonstrated via examples. 


TABLE I 

Comparison of computation times and operations eor (P, r)-MRPI and e-mRPI sets. 




LPs solved Minkowski 

sums 

Mean time (s) 

r = 6 


K = 

[-0.4345,-1.0285] 

1 

0 

0.005 

r = 20 



1 

0 

0.007 

r = 48 



1 

0 

0.019 

e-mRPI 01 (r = 

48) 


369 

11 

2.9 

r = 20 


K = 

[-0.0796, -0.4068] 

1 

0 

0.008 

r = 60 



1 

0 

0.036 

r = 172 



1 

0 

0.30 

e-mRPI 01 (f = 

172) 


3250 

42 

25 
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Fig. 2. Comparison of (P, r)-mRPI and e-mRPI sets for the system l ll5t with K = [—0.4345, —1.0285] and different disturbance 
sets. 
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